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Skew orthogonal polynomials arise in the calculation of the n-point distribution function for 
the eigenvalues of ensembles of random matrices with orthogonal or symplectic symmetry. 
In particular, the distribution functions are completely determined by a certain sum involv- 
ing the skew orthogonal polynomials. In the cases that the eigenvalue probability density 
function involves a classical weight function, explicit formulas for the skew orthogonal poly- 
nomials are given in terms of related orthogonal polynomials, and the structure is used to 
give a closed form expression for the sum. This theory treates all classical cases on an equal 
footing, giving formulas applicable at once to the Hermite, Laguerre and Jacobi cases. 
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1 Introduction 

The classical polynomials play an essential role in calculating statistical properties of the eigen- 
values of certain classes of random matrices with normally distributed elements. As a concrete 
example, consider a random N x N Hermitian matrix X in which the diagonal elements (which 
^ . must be real) and the upper triangular elements are independently chosen with normal distri- 

butions N[0, l/\/2] and N[0, 1/2] +zN[0, 1/2] respectively. Then the corresponding eigenvalue 
probability density function (p.d.f.) is proportional to 

N 

Y[w 2 (x j ) Yl (xk-Xj) 2 , w 2 {x) = e- x \ (1.1) 

i=l ^<j<k<N 

2 

Note that this involves the classical weight function e~ x . 

The monic orthogonal polynomials corresponding to the weight function W2(x), {pj(x)} say, 



are introduced by writing the product of differences in (1.1) as the square of the Vandermonde 
determinant, and adding together appropriate multiples of the columns. This shows that (1.1) 
is equal to 

N 



j,k=l,...,N 



2 

(1.2) 



This form of the p.d.f. provides the starting point for the calculation of the n-point distribution 
functions (see Section 3), which in turn are used in calculating spacing distributions (see e.g. [||). 

The p.d.f. fll.lD with the classical weight function W2(x) = x a e~ x , (x > 0), and 102(2;) = 
x a (l — x) b , (0 < x < 1), also occurs in random matrix problems (see e.g. ||). Thus if A is a 
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rectangular n x m, (n > m), matrix with complex elements independently distributed according 
to the normal distribution N[0, l/\/2] + iN[0, l/y/2], then the eigenvalues of X^X have the 
distribution with ii^ (») = x n ~ m e~ x and N = m. Similarly, if X\ and A2 are rectangular 
matrices of dimensions n% x m and ri2 x m, then the eigenvalues of x|Xi(l + x\X2)~ l have the 
distribution ( |P] ) with W2(x) = x ni_m (l — x) n2 ~ m and N = m. 

If the random matrices X which lead to the eigenvalue p.d.f ( |1,1[ ) with the classical weight 
functions noted above are specified to have real elements, or real quaternion elements represented 
by 2 x 2 blocks of the form 

(in which case the eigenvalues are doubly degenerate), the corresponding eigenvalue p.d.f. is 
then proportional to 

N 

Y[w^(xj) Y[ \x k -Xjf. (1.3) 

j=l l<j<k<N 

Here (3 = 1 when the elements are real, while (3 = 4 when the elements are real quaternion, and 
the weight function is again of the same classical form as in the complex case ( |1 . 1[ ) (up to scaling 
of the eigenvalues Xj). We remark that ensembles of such matrices have orthogonal {(3 = 1) and 
symplectic {(3 = 4) symmetry. Instead of writing the product of differences in (|1.3j ) in terms of 
an ordinary determinant involving skew orthogonal polynomials, it has been shown by Mahoux 
and Mehta § (see also ||, §, |l|]) that for the purposes of calculating the n-point distribution 
one should introduce a quaternion determinant involving skew orthogonal polynomials. The 
definition of a quaternion determinant can be found for example in Q|; we remark here that 
for the 2N x 2N matrices Q occuring in such a calculation (which technically are self dual) the 
quaternion determinant, to be denoted Tdet, is related to the usual Pfaffian via the formula 

TdetQ = Pf(ZQ), Z:=1 N ®( ~ 1 j . (1.4) 



Our interest is in the calculation of the classical skew orthogonal polynomials in terms of 
classical orthogonal polynomials, and the application of these formulas in the evaluation of 
the matrix elements in the Tdet expression for the n-point distribution. We begin in Section 
2 by revising the relationship between the p.d.f. ( |1.3[) for (3 = 1 and 4 and skew orthogonal 
polynomials. Also revised is the corresponding expression for the n-point distribution functions 
in terms of the skew orthogonal polynomials. Next we present the recent theory of Adler and 
van Moerbeke [Q] which identifies the family of orthogonal polynomials naturally related to a 
particular family of skew orthogonals, and furthermore relates the skew orthogonal polynomials 
at f3 = 1 to those at (3 = 4. 

In Section 3 this theory is used to calculate the classical skew orthogonal polynomials in 
terms of their natural basis of orthogonal polynomials. The results are shown to be consistent 
with expressions obtained in the earlier studies of Nagao and Wadati [12|. Application of these 
expansion formulas is made in Section 4, in which the fundamental series involving the skew 
orthogonal polynomials occuring in the expression of the n-point distribution is summed for 
all the classical ensembles. In the Hermite case this reclaims known results. For the Laguerre 
ensemble our summation formula differs from that obtained in recent studies of Widom 14] 
and Forrester et al. || — the reconciliation of the two seemingly different forms is made in the 
Appendix. The remaining classical case is the Jacobi ensemble. Here our summation formulas 
constitute new results. 
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2 Revision 



2.1 Skew orthogonal polynomials and the n-point distribution 



The n-point distribution pr n \ with respect to the p.d.f. ( |1.3j ) is defined as a ratio of multiple 
integrals according to 



P(n){Xl,. . . ,. 



3=1 



l<j<k<N \ x k — Xj\ 
l<j<k<N \ x k — 



(2.1) 



Generalizing the pioneering work of Dyson [Q], it has been shown by Mahoux and Mehta |?J that 
for [5 = 1 and (5 = 4 /j( n ) can be evaluated as an n x n quaternion determinant involving certain 
skew orthogonal polynomials. 

In preparation for presenting these results, we recall that an inner product (/, g) is referred 
to as skew if 

(f,g) = -(g,f) 

while a family of (monic) polynomials {q n (x)}n=o,i,2,... are said to be skew orthogonal if 

(q2m,Q2n+l} = ~ (</2n+l, 0.2m) = r m 5 m ,ni (<?2m, Q2n) = (<?2m+l > Q2n+1 ) = 0- (2-2) 

Equivalently, for any given iV we require 

n 
-n 



[(Qj>Qk)]j,k=0,...,2N- 



R, R '.- 



r 2 
-r 2 



r N 
-r N 



(2.3) 



We remark that the relations ( |2.2[ ) are unchanged under the replacement 

q2m+l{x) l-> g2m+l(^) + a 2 m<72m0z) (2.4) 

for arbitrary a2 m , so the skew orthogonal transformations are non- unique up to this mapping. 



The specific skew inner products of relevance to the calculation of (2.3) are 



and 



1 f°° 

</,<?>! := 2 j^e-VM I dye- v ^sgn(y-x)f(x)g(y) 



(f,9h ■= - I Jxe- 2V ^\f{x)g'{x) - f'(x)g(x)) 

i) we have chosen 



2 J — OO 

for (3 = 1 and /? = 4 respectively, where in (^ 

-V(a:) 



-2V(x) 



(2.5) 
(2.6) 

(2.7) 



We denote the corresponding monic skew orthogonal polynomials by {9^}j=o,i,... an d {^j 4 "^}i=Q,i,- 
and apply similar superscripts to the matrix R and its elements. 

With these preliminaries, the quaternion determinant formulas for ( |2.1| ) can now be pre- 
sented. 
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P — 1, N even 

For (3 = 1 the result depends on the parity of N . Let iV be even and set 

®k{x) := - / sgn(x - y)qi'(y)e- viy) dy, 

£ J —CO 

h(x,y) = 



with 



Si(x,y) h(x,y) 
Di(x,y) Si(y,x) 



N ^ e -V{y) , (1) 

Si(x,y) = 22 (1) ($2k(x)qi k > +1 (y) - $2k+i( x mk(y) 

k=o r k 
d 

Di{x,y) = —Si(x,y) 

If 00 1 
h[x,y) = -J S 1 (x,z)sgn(z - y)dz - -sgn(x - y) 

Then with w\{x) = e~ v ^ in ( |2~T| ) we have 

P(n)(xi, ...,X n )= qdet[/i(Xj,X fc )]j,fc=l,...,n- 



(2.8) 



(2.9) 
(2.10) 



/3=1, N odd 



Let {qn }n=0,l,... and {r^}^ 



,(jV-l)/2-l 



be as in the iV even case, and define 



fodd 
Jl 

Sf d (x,y) 

Df d (x,y) 
If d (x,y) 



\ roa 

2 



dxe-^q^ix) 



(N-i)/2 ■ 

q^(x) := (|_°° dx'e-^q^H^y^x) (n = 0, . . . , N - 2) 



1 f 00 

(iy e _y ^- ) sgn(x — (x) (n = 0, . . . , N — 1) 



2 7 — oo 

^ dd (x,y) I? dd (x,y) 
Df d (x,y) Sf d (y,x) 

1 v-m- 1 e -V(y) 



E 

k=0 



.(1) 



(*2fc(a:)5£ ) +1 (y) " ^2k+i{x)q^(y) 



(i) 

(jV-l)/2 



+ 



2r 



9a; 



^ odd (x,y) 



1 /"°° i 1 

- J S° dd (x, z)sgn(z - y)dz - -sgn(x - y) + 



$JV-l(s) 

2r (i) 



With this notation, ( 2. ICQ applies with /i replaced by /° dd H. 



(2.11) 



= 4 

For (3 = A the results are independent of the parity of N . Let 



iV-l e _y(x) 



in i) 2rm^ 



dy 



dy 



4 



d 

D 4 (x,y) = —S 4 {x,y) 



h(x,y) :-- 
and set 



S A (x,y')dy' 



(2.12) 



S 4 (x,y) h(x,y) 
D 4 (x,y) Si(y,x) 



h{x,y) :-- 

Then with w 4 (x) = e~ 2V(x) in (ED]) we have @, 

P(n)(xi, ...,x n ) = Tdet[/ 4 (zj,a; fe )] Jife=1) ... )r , 

2.2 Construction of the skew orthogonal polynomials 



(2.13) 



A number of works have addressed the construction of the skew orthogonal polynomials 0, 12, 
[|, ||. In particular Nagao and Wadati [jl2| have succeeded in giving explicit formulas for the 
classical skew orthogonal polynomials in terms of related orthogonal polynomials. More recently 
Adler and van Moerbeke 0J, generalizing some earlier work of Brezin and Neuberger Q, have 
identified a natural orthogonal polynomial basis for the skew orthogonal polynomials {q n (x)} 
corresponding to the skew symmetric inner products (|2.5| ) and (|2.6| ) and shown how to use this 
basis to specify the q n (x). Here we will revise this formalism, which in the next section will be 
implemented in the classical cases and the results compared with those known from the work of 
Nagao and Wadati Jl2[. 



Following @, with V{x) as in (EDI), let 



2V'{x) 



g(x) 

fix) 



be a rational function, and introduce the operator 

d , ff'-9 



n 



J dx V 2 ) 
Furthermore, introduce the symmetric inner product 



~ 2V ^(/)(x)ip(x)dx. 



(2.14) 



(2.15) 



(2.16) 



Then it is shown in 01 that 



!>,n _ V)2 



-(n(t),ip) 2 , 



V(x)^V{x)+\ogf{x)' 
V(x)^V(x)-lo S f(x) 



(2.17) 



(the first of these results was noted earlier in [Q]). 

The practicality of these formulas lies in the fact that with {pj(x)}j = o,i,... denoting the 
monic orthogonal polynomials corresponding to the symmetric inner product ( 2. 16] ), the matrix 
[(pj, n>Pk)2]j,k=o....,N-i (which by the first formula in ( 2.17 ) is skew symmetric) has non-zero 
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terms only in a finite number of diagonals about the main diagonal for / and g polynomials in 
(2.14). In particular, for the classical polynomials this number is precisely two, giving 



M := [(pj,np k ) 2 ]j,k= 



1,...,N-1 



' 


co 








-co 





ci 








-ci 





C2 








-C2 






c N -i 
-cjv-i 



(2.18) 



Equivalently, for {pj(x)}j=o,i,... a se t of classical orthogonal polynomials, 



np k (x) 



Ck 



(Pk+i,Pk+ih 



Pk+i(x) + 



Cfc-1 



(Pk-l,pk-l) 



-Pk-i{x). 



(2.19) 



The simple structure of ( |2.18j ) allows the classical skew orthogonal polynomials to be ex- 
pressed in terms of their orthogonal polynomial counterparts. Let {qj (x)} and {qj 4 \x)} 



denote the skew orthogonal polynomials corresponding to the inner products ( |2.5D and (]2 
with V(x) i— > V{x) + log/(x) and V(x) t— > V(x) — log /(a;) respectively. Similarly let {rj 1 ^} 
and {^j 4 ^} denote the corresponding normalizations as in fl2.3|) . We can write the matrix 



i=o,...,iV-i 



[(pj,n 1 Pk)2]j,k=o,...,N-i in terms of M. Thus with 

D := diag((p i ,p i ) 2 
we have 

N-l / _i \ 

([(p j ,n _ V)2]j,fc=o,...,Af-i-D _1 AA) = - V 2 (pi,np k ) 2 

where use has been make of the first equation in ( 2.17| ) . Now according to ( 2.1S| ) 

N-l i \ 

E {P ( UnPk ^ Pl{x)=np k {x) for k^N-1. 

Substituting this back in ( |2.2l| ) shows 

[{pj,n~ 1 p k ) 2 ]j i k=o,...,N-iD' 1 N' = D' 



(2.20) 



(2.21) 



(2.22) 



where D' = D + C k with C k having all entries zero except for the final column (the explicit 
value of these entries will not be required) . 

Next introduce transition matrices QW and such that 



[</,• (x)] j=0 ,...,N-i = Q^\Pj(x)] j=0 ,...,N-l, 



s = 1 or 4. 



(2.23) 



Both these transition matrices must be lower triangular matrices with l's along the diagonal. 
Then we have 



V(x)^V(x)+logf(x) 



Q (1) [(Pj,Pk)l v{x) ^ vix)+logf{x y 

-Q( 1 )[( Pj ,n-V)2b, fe =o,...,iV-iQ {1)T = -QVd'N^DQW (2.24) 



j,k=0,...,N-l 
]j,k=0,...,N-lQ^ T 
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where in obtaining the first equality of the final line the second formula of ( [2.17D has been used, 
while in the final equality (2.22) has been used. Taking inverses of both sides of ( |2.24| ) shows 
that it is equivalent to the statement that 



Q^ T R^- l Q^ = -D^ND 



1-1 



A similar calculation shows 



-(4)> 
Wj ,Q k /4 



V(x)^V(x)-\ogf(x) 



j,k=0,...,N-l 



or equivalently 



(2.25) 

(2.26) 
(2.27) 



3 The classical skew orthogonal polynomials 

3.1 Summary of known results 
Hermite case 

Here e~ v ^ = e~ x2 ' 2 . The monic orthogonal polynomials with respect to the inner product 
fl2.16j ), and corresponding normalization are then 



Pk\x) 



2- k H k (x) 



(3.1) 



The skew orthogonal polynomials and corresponding normalizations can be expressed in terms 
of these orthogonal polynomials according to 0] 



and 



-(i) 



l2m+l( X ) 
4m( x ) 



P2m(x) 

P2m+l(x) - mp 2 m-l(x) 

2- 2n V^T{2n + 1) 



P2m+l(x) 
m 

ml —,P2n(x) 



n=0 
-2n-l 



V7rr(2n + 2). 



(3.2) 



(3.3) 



Laguerre case 

Here we have 



e -V(x) 
Pk(x) 

(j>k,Pkh 



x a/2 e -x/2 



(-l) k k!L%(x) 

r(fc + l)T(a + k + 1) 



where Lf. denotes the Laguerre polynomial. With 



(3.4) 



-V(x) 



-V(x)-logf(x) 



x (a-l)/2 e -x/2 



(3.5) 
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in ( |2,5D the corresponding monic skew orthogonal polynomials and normalization are [jl^ l 

& ( X ) = P2m(x) 

$bi(x) = -(2m + l)!L^ 1 +1 (x)-(2m)!(a + 2m)^ m ( 3 ;) 

= P2m+i (x) + (2m + l)p2m (x) - 2m(a + 2m)p 2m -i (x) 
f$ = 2r(n + l)r(a + 2n + 2) (3.6) 



where in obtaining the second equality in (|3.6| ) use has been made of the formulas 

L a n - 1 (x) = L a n (x)-L a n _ 1 (x) ±L%{x) = -L«±\(x). (3.7) 
Recalling the fact that the skew orthogonal polynomials are non-unique up to the transformation 



(2A) we see that we can equally as well write in place of the second formula in (J3 



<?2m+i( x ) = P2m+i{x) - 2m(a + 2m)p 2m -i{x). (3.8) 

At = 4 put 

e -W(x) _> ^(V^-log/^)) = x a+l e -x (3 9) 

in fl2,6|). Then we have O] 

4 4 2 +1 (x) = p 2m+ i(.) 

m -i 

£(*) = 2 2 -m!(a/2 + m)!^ 22jj!(Q/2+j)! ^(x) 

fW = ir(2n + 2)r(a + 2n + 2). (3.10) 

Jacobi case 

In the Jacobi case 

e - v ^ = (l-x) a/2 (l + x) b/2 . (3.11) 
The monic orthogonal polynomials and normalizations with respect to ( 2.16| ) are then 



v,(x) - ofcM r(a + b + fc + l) (g , 6) 

, s o(a +fe+ i+ 2fc ) r(fc + i)r(a + 6 + fc + i)r(a + 1 + k)r(b + i + fe) 

(M)2 = 2 r(a + 6 + 2fc + 2)r(a + 6 + 2A: + l) (3 " 12) 

With 

e -V(x) _> e -V»-log/(*) = (1 _ x )(a-l)/2 (1 + x) (b-l)/2 (3 _ 13) 

in ( |2.5[ ) the skew orthogonal polynomials and corresponding normalization at p = 1 are |jX 

= ^(2m)I ^ t 6 1 4m + lj { (2m + 1)(Q + b + 2m " ^t^"^) 
(a + b- 2)(q + 2m)(6 + 2m) d ( -i,6-i), 
(a + 6 + 4m)(a + 6 + 4m + 2)dx 2m 

= p 2m+1 (x) + 2(2m + 1) — + Q ~ + 6 b - 2 p 2m (x) 
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(2m + a)(2m + 6)(a + 6 + 2m) 
~ m (a + 6 + 4m - l)(a + 6 + 4m) (a + 6 + 4m + l)(a + 6 + 4m + 2 ) P2m ~ 1 ^ 
~(i) _ oa+fe+4n+2^, r(a + 2n + 1)T(6 + 2n + l)r(a + 6 + 2n + 1) 
" " [ r(a + 6 + 4n + l)r(a + 6 + 4n + 3) 1 j 

where in obtaining the second equality in the formula for q2m+l{ x ) use nas been made of the 
formulas 

(2n + a + + l)P^\x) = 

(n + a + + + a + /3 + 2) (a+1|/?+1) 

( 2n + a + + 2) ^ n {X> + 

(n + a + /? + l)(n + /3+ 1) {n + a + + \){n + a) \ {a+lfi+l) 
{2n + a + + 2) + {2n + a + 0) ) n ~ x W 
(n + a)(n + 0) » 
(2n + a + /3) ™" 2 W " 

The non-uniqueness of the skew orthogonal polynomials up to the transformation (|2.4|) implies 



we can equally as well write q~2m+i(%) in ( 3.14 ) as 

<72m+lO) = P2m+l( a 

_ 8 (2m + a)(2m + 6)(q + 6 + 2m) 

(a + 6 + 4m- l)(a + 6 + 4m)(a + 6 + 4m + l)(a + 6 + 4m + 2) P2m_ll2;j ' 1 j 

At /3 = 4, make the replacement 

e -2V(*) _> e -2(V( X )-lo S f( X )) = {1 _ x) a+l {l + x) 6+l ( g lg) 

in fl2.6| ). The corresponding skew orthogonal polynomials and normalizations are then 



4m( X ) 



2 6 m+ a+b m iY( a /2 + b/2 + m + l)r(o/2 + m + l)r(6/2 + m + 1) 
? r(a + fe + 4m + 2) 

v 1 r((a/2 + fe/2 + j + l) r(a + fe + 4j + 2) 
x ^ r(o/2 + j + l)r(6/2 + j + 1) T(a + 6 + 2j + l) P2j (xj 

2 6 m r(o/2 + 6/2 + m + l)r(a/2 + m + l)r(6/2 + m + 1) 
m ' T(a + 6 + 4m + 2) 

^J- r(a + fe + 4j + 2) 

X j\2V r(o/2 + 6/2 + j + l)r(a/2 + j + l)r(6/2 + j + if 2 ' W 

(4) _ 2 a+fe+4 "+ 2 r(2n + 2)r(a + 2n + 2)r(6 + 2n + 2)T(a + 6 + 2n + 2) 
r " ~ T(a + 6 + 4n + 2)r(a + 6 + 4n + 4) ^ ' ' 

where in obtaining the second equality in the formula for q^l_(x) use has been made of the 
duplication formula 

T(z)T(z + 1/2) = 2 1 / 2 - 22 (27r) 1 / 2 r(2z). 
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3.2 Derivation from general formalism 

= 1 



The expression for skew orthogonal polynomials at (3 = 1 is determined by solving the equation 
(PI ) for Q« = [Q jk ] jtk =o,...,N-i- In fact this can be done by premultiplying both sides by 
(Q^ T ) _1 , and equating the lower triangular entries on both sides (not including the diagonal). 
Using (p. 3[) we see that on the l.h.s. 



(1) 



Q 2 oM 1) 
-Qm/'> 



(i) 

2 



-Q3iM 1] 

-Q5i/rf } 

Q4l/^ 1} 



1/ff 
Q42/^ 1} 



1/f 



V 



/ 



2iVx2iV 

(3.18) 

where the subscript in ()_ denotes the strictly lower triangular entries. On the r.h.s. we note 
from the fact that (Q^ T )~ 1 is an upper triangular matrix with l's along the diagonal and the 
explicit form ( gig ) of M that 



((Q« r )- 1 irW 1 



( * 
7o 



V : 



71 




72 



(3.19) 



where 



7j : = Cj/{p j+ l,p 1+ l) 2 {Pj,Pj)2- 



Equating ( 3.1§| ) and ( |3.19| ) gives 



72p 
Q2p,l 
Q2p+l,l 
Q2p+l,2p-l 



o (z = 
o (/ = 

-72 P -r 



(p = 0,...,JV 
0,...,2p-l) 
0,...,2p-2) 



1) 



-72p-i/7: 



2p 



(3.20) 



(3.21) 



while Q2p+i,2p is left unspecified. Hence 



P 



P2 



^5+1 ( ffi ) = P2j+\( X ) + Q2p+l,2pP2j {X 



72j 



■l/72p. 



(3.22) 



Note that the fact that C?2p+l,2p is arbitrary in the formula for q^ +l {x) is consistent with the 



skew orthogonal polynomials being non-unique up to the transformation (2.4). The simplest 
choice is Q2 P +i,2p = 0. 
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Let us verify that fl3.22|) with Q2 P +i,2p = reclaims the results . 1|) , (|3.6D and (|3,14|) . First 
7j must be specified (recall ( |3.20| )). From fl2.19| ) and ( |3.20 ) we have that 



nPk(x) = -7k(j>k,Pk)2Pk+i(x) + lower degree term, 



(3.23) 



which together with the explicit form of n for each of the classical weight functions and the fact 
that each pk{x) is monic implies 



lk(Pk,Pk) 



1, 
l 

2" 



Hermite 
Laguerre 



(3.24) 



\{2k + 2 + a + b), Jacobi 



Furthermore, the explicit form of (pk,Pk)2 is given by ( |3.1| ), ( |3,4j ) and (3.12) in the Hermite, 
Laguerre and Jacobi cases respectively. A straightforward calculation then enables the formulas 
of the previous section for f3 = 1 to be reclaimed from the general formulas ( |3,22j ). 



P 



At P = 4 we must solve the equation ( 2.27j ) . This we will do for the matrix (Q^ 4 )) -1 := L = 
lPjk]j,k=o,...,N-i, where L is lower triangular with l's along the diagonal. In this case it is 
sufficient to multiply both sides of ( |2 . 2 7 ) by Q^ T and equate the strictly lower triangular parts 
of both sides. On the l.h.s. we have 



/ 



\ 



f (4) 
'0 P21 

P31 



ft 



41 



'0 P20 

r l P30 
'1 P40 



_ f (4) 

=(4) 1 3(4) ~(4) S(4) 



'0 



43 



ft 



12 



(3.25) 



2Nx2N 



while on the r.h.s. 



(NQ 



(4)T 



(JV). 



(3.26) 



Equating ( |3.25| ) and ( |3.26j ) gives 

C2p = 

/3 (4) 

P2p,2p-2 — 



f(4) 
P 






-Pap-i/r^Zi 



p = 0,...,iV-l 
j = 0,...,2p-l 
j = 0, . . . , 2p - 3, j 

-C2p-l/c2p-2 



2p- 1 



(3.27) 



while 2p i s undetermined. Thus we have 



P2j+l(x) 
P2j(x) 



it^+ftit^ix) 



f(4) 
'p 



if(x) 



C2p- 



C 2 j-2 



(3.28) 
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lm- 



The non-uniqueness of the skew orthogonal polynomials up to the transformation (| 
plies we can choose /?2j+i 2j = ®- Doing this, and solving for q£j (x) in the second equation of 
( |3!2SD gives 



f(4) 
P 



P2j+l( x ) 

/i c 2p ^to K p=o C2 p 

C2p- 



P2l(x) 



(3.29) 



We recall that the quantities Cj are, according to ( |3.20| ) and ( |3.24 ), simply related to the norm 
(Pj+iiPj+i)2 which in turn is given by ([0]), ( |3.4| ) and ( 3.12 ) in the Hermite, Laguerre and Jacobi 
cases respectively. Prom this it is straightforward to verify that the general formula reproduces 
the results of the previous section for = 4. 



4 Summation formulas 

The results of Section ^O] show that the n-point distributions are given in terms of the quantities 
Si(x,y) and Si(x,y) (recall ( |2.9j ) and ( 2.12; )) at (3 = 1 and (3 = 4 respectively. It is of interest 



to compare these formulas with the expression for the n-point distribution (2.1) at /3 = 2. With 
W2(x) = e~ 2V ^ it is an easy consequence of (1.2) that 



P(n) (xi,...,x n ) = det[S 2 (x j ,x k )] j>k=lt ^ n , S 2 (x,y) := ^ T^T ■ C 4 - 1 ) 

Thus again the n-point function is completely determined by a summation. In fact for {Pj(a0}_?=o,i, 
a general set of monic orthogonal polynomials this sum can be evaluated exactly according to 
the Christoffel-Darboux formula 

pi(x)pi(y) = i PN(x)pN-i{y) - pn-i(x)pn(v) , 42 x 

^ (pi,Plh (p N -i,PN-i)2 x-y 
Here we seek an analogous summation formula for Si(x,y) and S4(x,y). 



There is some recent literature on this problem. Using different methods, Widom [14] and 
Forrester et al. || have expressed S\(x,y) and S^{x,y) in terms of the Christoffel-Darboux 
summation fl4,2j ) in which {pj(x)} are the monic orthogonal polynomials associated with the 
weight function W2(x) = e~ 2V ^ x \ plus a 'correction' term. In general the correction term (which 
is given in a different form in the two different formalisms j[4|] and ||) does not appear to 
have a closed form evaluation. Exceptions are the Hermite (H) and Laguerre (L) cases. Then 
the correction term factorizes as a function of x times a function of y, both of which can be 
computed exactly (the results in the Hermite case were known from earlier work; see e.g. [||). 
Explicitly 

S[ c \x,y) = Sf\x, y) + 2iV+1 ^ ( ^ _ iy H N ^(y) £^ sgn(x - t)e~ t2 ' 2 H N {t) dt 
Sf^\x,y) = s[ G \x,y) + Htf ~^ 



JZoe~ t2/2 H N -i(t)dt 
2Sf\x,y) = sP(x,y)\ N ^ 2N + e-y 2 / 2 H2 N (y)^= f ^J_^e- t2 / 2 H2 N - l (t)dt 
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s[ L \x,y) = S ( 2 L \x,y) + 
d 



-y a l 2 e-y/ 2 



2Sf\x,y) 



4T(N + a) 

{^ L N(y)) «gn(x - y)y a / 2 - 1 e-y/ 2 (L a N (y) ~ ^_i(y) ) dy 
S ( 2 L \x,y) 



+ 



N^2N 

(2N)\y a ' 2 e-yl 2 L% N (y) - L$ N _ x (y) f* . 



2T(2N + a) 



y 



d 



t a/2 e -t/2 L a 2N {t)dt. 

at 



(4.3) 



One feature of these formulas is that the quantities S±(x,y) and S<i(x,y), in which the skew 
orthogonal polynomials are with respect to the skew inner products (|2.5| ) and ( |2.6| ) (and thus 
involve e~ v ( x ' and e~ 2V ^ respectively) are expressed in terms of the weight function e~ 2V ( x ' . 
However, we have seen in the previous section that it is the skew orthogonal polynomials at 
(3=1 with the replacement 

V{x)^V{x)+\ogf{x) (4.4) 



m 



m 



5|), and the skew orthogonal polynomials at f3 = 4 with the replacement 

V(x) i ^ V(x) — logf(x) 



(4.5) 



.6|) which are naturally expressed in terms of the orthogonal polynomials for the weight 
function e~ 2V( - x \ This suggests that we consider the quantities Si(x,y) and S±{x,y), defined 
as in Q2.9| ) and ( |2.12|) but with the corresponding skew orthogonal polynomials defined with 
the replacements fl4.4|) and ( |4.5| ) respectively, and seek to sum them in terms of the orthogonal 
polynomials for the weight function e~ 2V<yX \ In the Hermite case nothing changes because, 
according to the definition ( 2.14| ), f(x) = 1 and so log f(x) = in (4.4) and (4.5). However, in the 
Laguerre case this is already a different viewpoint as then V(x) is modified by a t—. ► a± 1. Indeed 
we find that the simple structures ( |3.22| ) and ( 3.28| ) giving the skew orthogonal polynomials for 



the modified V(x) in terms of the orthogonal polynomials for the weight function e 2V ( X ) allows 
S±(x,y) and S^{x,y) to be summed in the classical cases by the same general formula in each 



case. 



(3 — 1, N even 

Let us write ^ 

Pi(x) = J2^^\x), #M = 1 (4.6) 
3=0 

where {gj (a;)}j=o,i,... is a set of monic skew orthogonal polynomials with respect to the inner 
product ( |2.5| ) (modified so that V is replaced by V), and {pz( x )}/=o,i,... is the set of monic 



polynomials with respect to the inner product ( 2.16 ). The quantity S\(x,y) can be expressed 
in terms of the polynomials pi(x). 



Proposition 4.1 Let N be even. With S\(x,y) defined by ft2.Q) , modified so that V i— ► V and 
q^\x) i— > q^\x), and {pj(x)}j=o i the set of monic orthogonal polynomials associated with 
the weight function e~ 2V( - x ^ (assumed complete) we have 

S 1 (x, s) = e-^H^)-n»)(^ 1 ^W + £ £ ^ q f\y)\ (4.7) 

\ n=0 [Pn,Pn)2 n=N k=Q {Pn,Pn)2 J 
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Proof Our derivation is motivated by the derivation of a similar formula in [11]. Now, from 
the definitions (p.8| ) and ( |2,5[ ) (the latter with V i— > V) it is possible to write 

$k(x) = e*&{$\y)\6(x-y)) 1 . 
But from the completeness of {pj{x)} we can make the expansion 

Pn(x)p n {y) 



[Pn,Pn)2 



Substituting 



for Pn{y) and using the skew orthogonality of {<?fc(y)} shows 

oo 



Pv\ X ) 5(1) 

>V2fc+l 



(j>v,Pu) 



2fc+l 



(*) 



v=2fc+l 

oo / \ 



x) f v~/ 



(4.8) 

(4.9) 
(4.10) 



u=2k 



Next substitute these formulas in the definition (2J)). The stated formula follows after minor 
manipulation involving further use of (4.C). 

The first sum in (4.7) is evaluated according to the Christoffel-Darboux formula (4.2). This 
applies for general monic orthogonal polynomials {pj(x)}. To evaluate the second sum requires 
knowledge of the transition coefficients . In the case of the classical polynomials these coef- 
ficients can be determined from ( 3.22[ ). Setting Q 2 p+l,2p = and comparing with the definition 
(4J3) we see that 



{1) 

P 2l,j 
^2l+l,2j 



s(i) _ nL=i a k 

P2l+l,2j+l 



0, (i = 0,...,2/-l) 
0, (j = 0,...,l) 

72fe-l 



Uk=l a k 

(1) 



a k :-- 



12k 



In fact we don't require the explicit form of the 0,- , but rather their factorization property 



50) - qO) 5(1) n> N 

Pnk — Pn,N-lPN-l.ki n — iV ' 



(4.11) 



which is evident from the above formulas. Substituting ( 4.11 ) in the double summation in (4.7) 
and recalling (|4.6[) shows 



oo N—l ( \ 

1^ 2s ~ \„ f 3 nk ( ik{y) 



n=N k- 



n 2 



Pn(x) 5(1) 
„ v (Pn,Pn)2 



Pn,N-l)PN-l(y) 



AD 

' N/2-l 



-V(*)+2V(*)$ N _ 2 ( x) 



PN-l(x) 
(PN-l,PN-lh 



PN-i(y) 



where the second equality follows from ( |4.9[) and the fact that N is assumed even. Thus 
Si(x,y) 

+ 7Ar_ 2 e" y(y) 4>Ar_ 2 (x)pAr_i(y) 



e -(V( X )-V(x)) e (V(v)-V(y)) S2 ^ y) 

e -<y{*)-v(*)) e <yiv)-v{v))s 2 ( x ,y) 



N^N-1 



N^N-1 



+jN-2e ^ (s/ W-i(y)^ 



u{x-t)p N - 2 {t)e' v{t) dt, 



(4.12) 
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where ")n-2 is specified in terms of (pn-2,Pn-2)2 by (|3.24| ). Hence the summation of Si takes 
the same form for all the classical cases. 

It is of interest to compare ( |4.12| ) with the summation formulas for Si in ( f4.3| ) . In fact neither 
the Hermite nor Laguerre summations in ( [4.3[ ) are of the same form as ( |4. 12|) . B ut in both cases 
the equality of the different forms can be established directly. Consider ( |4. 12 ) in the Hermite 
case when it reads 



s{ G) (x,y) = sf\x,y) 



ivwiv-i 1 2 n tt 1 / 2 (N - 2)! 



+ 



dt sgn(x - t)e~ t2/2 H N ^ 2 {t). 



This can be brought into agreement with the formula for §i in ( |4.3| ) if we note from the 
identity 



H 2k (y) 



l 



that 



$2k(x) 



4(ft + l/2) 
1 



H- 



2k+2 



(y) 



-2k 



-<Z>2k+2(x) 



-e- x */ 2 H 2k+ i(x). 



(4.13) 



k + l/2" K ' r ' K ' 2k + 1 

In the Laguerre case the manipulation of (4.12) into the form given in ( |4.3j ) is undertaken in the 
Appendix. 



13 = 1, N odd 

In the case of N odd, we see from the definitions ( 2.11| ) that 



Sf d (x,y) = Si(x,y) 



+ 



£li(y) 



+ 



N^N-l 2sAf_i 



Yl ..Hi 



SN-l 



s2kQ2k+i(y) + s 2 k+ihk{y) 



SN-l 



E 6 Jl) ( ~ S2k®2k+l(x) + S2k+l®2k(x, 



k=o r k 



where 



Sk := 5 



-v^$\ x )dx. 



(4.14) 



The quantity Si(x,y) is evaluated by ( [4.12 ). Furthermore, from the definitions we see 

that 



(iV-l)/2-l v(y) 
E 



fc=0 



S2kq2k+i{y) + s 2 k+iq2k(y) 



lim Si(x,y) 



N^N-l 



(4.15) 



where here the second equality follows from (4.12), while 



(AT-l)/2-l 



1 



E 7(1) ( ~ Slk$2k+l(x) + S 2 k+l®2k(x i 
k=0 f k 



1 f°° ( ~ 

- / sgn(x-y) - lim Si(x',y) 



N^N-l 



dy = ~'jN-3S N - 3 (f> N -2(x) (4.16) 
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where here the second equality follows from ( 4.15 ) and 

1 



4>j(x) :-- 



Hence 



^ dd (x,y) = S 1 (x,y) 



+ 



v ^sgn(x - y)pj(y) dy. 



PN-i(y) 



(4.17) 



N^N-l 2s N -i 



-V(y) 



-"/N-3SN-3- 



Y((f>N-l(x)pN~2(y) ~ PN-l(y)4>N-2(xfj , (4.18) 



where the fact that q$_i(y) = pN-i{y), which follows from ( 3.22 ) and the fact that N is odd, 
has been used. 



Let us show how to write ( 4.1Sj ) in the Hermite case as presented in ( |4 . 3| ) . Now, from the 
definite integral 

e~ x /2 H n (x) dx = (27r) 1/2 2" /2 (n - l)(n - 3) • • • 3 • 1 n even 



we see that 



SN-3 n 7AT-3 
lN-3- = 2 — ^ = 7iV-2- 



SN-1 



N -2 



(4.19) 



This fact shows that the final term in the second line of ( 4.18 ) is equal to the final term in the 
formula ( 4.12| ) for Si(x,y). Also, making use of the first equality in ( 4.1 9| ) and ( (4,13| ) shows 



~V^p N _ 2 (y)~f N _ 3 (^) N _ 3 (x) - ^-JL(f> N ^i(x) 



SN-1 



-V(x)-V(y) PN-2(x)pN-2{y) 



{PN-2,PN- 



2 2 



We see from ( 4.1 2| ) with N ^ N — 1 that this identity provides the final step in identifying 
fl4.18j) in the Hermite case with the formula in (4.3). 



= 4 

The strategy for obtaining a summation formula at (3 = 4 is similar to that used at (3 = 1. We 
write 

Mz) = E/fte, = l (4-20) 

3=0 

where {<jj^ (%)}j=o,i,... is a se t of monic skew orthogonal polynomials with respect to the inner 
product ( pOI ) (modified so that V is replaced by V) and {pz(^)};=o,i,... is the set of monic 
polynomials with respect to the inner product ( |2.16| ). 

Proposition 4.2 With Si{x,y) defined by modified so that V 1— ► V and q£\x) \— > 

q^\x), and {pj(^)}j=o,i,... the set of monic orthogonal polynomials associated with the weight 
function e~ 2V ^ (assumed complete) we have 

1 - - /2JV-1 / n / s 00 2JV-1 / x N 

S i {x,y) = \e-^-v(*)+ny)( £ £ £ -/^>-&q k (x) ) . (4.21) 

2 V ^ (Pn,Pn) 2 n f^v ^0 (^iV>2 J 
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Proof Again our derivation is motivated by the workings in [ 11 1 . Prom the fact that 



,V(x) 



6 dx Q2m[X) + dxV 92m[ 



d 



X 



we can check from the definition ( |2.6j ) that it's possible to write 

ii^&i-)) = ^ (a ° - y)\&(y)U- 

Substituting ( |4.8| ) for 5(x — y), and then substituting (4.20) for Pj(y) and making use of the 
skew orthogonality of {q^ix)} shows 



d 
dx 



Pu(x) s(4) 



u=2m+ 



(4.22) 



Proceeding similarly we can also show 



d_ 

d:X 



OO / \ 

fp-^W,j( 4 ) ( T \\ - f (4) p V(x)-2V(x) \p Pv\ x ) 5(4) 



u=2m 



(Pu,P, 



v)2 



(4.23) 



Apart from the sign, these equations are formally the same as those in (4J}) and ( 4.10 ). The 
stated formula thus follows as in the derivation of ( |4.7| ). 

Since the first sum is evaluated by the Christoffel-Darboux formula (4.2), it remains to 
evaluate the second sum in ( [4.21|) . This in turn requires the value of the p^. Now, according 



to ( 3.28 ) with /32p+i,2p = 0> in the classical cases the only non-zero value of (3^ for n > k is 



(-1) 



/ J 2n,2n-2 



C2n-1 
C2n 



where c n is specified by ( |3.20| ) and ( 3.24 ). Hence in the classical cases 



S 4 ( x , y ) = L e -(v(y)-v(y)) e v(x)-v(x) S2{x7y) 



N^2N 



1 

+-e 



-V(x)+V(y)-2V(y) c 2N-l P2N(y) (4) , , 
I \ 12N~2\- L )- 

C2N (P2N,P2N)2 



(4.24) 



We note that reference to the skew-orthogonal polynomial q^_ 2 ( x ) can eliminated by noting 



from ( |4.22| ), and the fact that the only non-zero value of P^l m +i is /^2m+i 2m+i = x ' that 



5(4) 



- V{x) q 2N -2(x) 



s( 4 ) 



r 



7V-1 



(P2iV-l,P2Af-l)2 A 



(4.25) 



Using ( [3.28 ) to eliminate fff_ l7 and recalling ( 3.20| ), from this formula we thus have 



Si{x,y) 



-(V(y)~V(y)) e V(x)-V(x) g , x 



N^2N 



\e^- 2V ^h2N-lP2N{y) r<H<-* ( "/) 2 .V-l(/). 



(4.26) 
(4.27) 



In the Hermite ca se (|4.26 ) gives immediate agreement with the formula in ( |4.3| ). However 
in the Laguerre case ( [4.26] ) has a different structure to the formula in ( f4.3| ). The manipulations 
necessary to show that the two formulas do indeed agree are undertaken in the Appendix. 
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5 Summary 



Classical skew orthogonal polynomials occur in the calculation of the n-point distribution func- 
tion pr n \ associated with the eigenvalue p.d.f. ( |1.3| ) with wp(x) a classical weight function and 
(3=1 and 4. In particular they occur in a certain sum, denoted Sp{x, y), which determines P( n ) . 
This sum is specified in fl2.9|) for /3 = 1 and iV even, in ( |2.11| ) for (3 = 1 and N odd and in fl2.12| ) 
for (3 = 4. The main achievement of this paper has been the closed form evaluation of Sp(x,y) 
in terms of particular classical orthogonal polynomials naturally related to the classical skew 
orthogonal polynomials. A self contained presentation of the relationship between the classical 
orthogonal and skew orthogonal polynomials is undertaken in Section 3. The summation for- 
mulas for S/3(x,y) (the tilde on Sp denotes that the weight functions are modified according to 
(4.4) and (|4.5| )), which apply equally as well to all the classical cases are given by fl4.12| ) for 
f3 = l and N even, (p^ ) for (3 = 1 and N odd and fl4~26p for f3 = 4. 
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Appendix 

Here the summation formulas of (4.3) for S/3{x,y) in the Laguerre case will be shown to give 
agreement with expressions implied by the general formulas ( [4.12 ) with the substitutions ( |3,24 ) 
and (3.4). Consider first the case (3 = 1. Because Si(x,y) is defined with the replacement (4.4) 
and f(x) = x in the Laguerre case, we see that we are required to show 

S 1 (x,y) = S 1 (x,y) (A.l) 

ay—>a+l 

We see from ( |4. 12| j and ( |3.4j ) that the r.h.s. is presented as a series in the linearly independent 
set of functions {y a ^ 2 e~ y / 2 (y)} . On the other hand, the l.h.s. can also be expressed in terms 
of this basis by using the first identity in ( |3.7| ) (with awa + l). Thus we must show that the 
coefficients of each y a ^ 2 e~ y ^ 2 L^ +1 (y) (which are functions of x) agree. The largest value of k 
which occurs is k = N — 1. On the l.h.s. the coefficient is then 



~— —Ln-i{ x ) - 7F7 ' An / sgn(x -u^^e-^L^i^du (A.2) 



.*,2„- X ,2 (N-l)\ N\ 

T(a + N) LN - l{x) -±T(a^N) ./„ 

while on the r.h.s. the coefficient is given by 

(N-iy i\ „,,.,,-, 

4r(a + N) Jo 
Now, by making use of the formula 



poo 

/ sgn(x - u)L a +\(u)u a / 2 e- u l 2 du. (A.3) 



I r((n+3)/2))r(q+n+2) 

a l 2 e- X ' 2 L« +l {x)dx= \ 2«/2-ir(n+2)r((n + a+3)/2 ) )> n eVen ' (A 4) 

0, n odd 
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we can check that both (|A.2| ) and ( A. 3 ) agree in the limit x — » oo. It thus suffices to show that 
both expressions have the same derivative. This we can do by making use of the identities (3.7) 
as well as the further identity 



xL a N + \(x) = -NL a +\x) + (2N + a)Lp-\(x) +-(N + a)L a ^ 2 (x). 



a+l 



0+1 



-a+l 



(A.5) 



The coefficient of y a / 2 e-y/ 2 L a k +1 (y) for k < N on the l.h.s. of flAp 

is given by 

m , (k + iy. 



T(a + l + k) 
while on the r.h.s. it is given by 



L%{x) 



k\ 



T(a + 2 + k) 



-a+l i 



T(a + 2 + k) 



xLj. ' ' {■)■) | -<■ 



a/2 -a/2 



(A.6) 



(A.7) 



Use of ( A.5 ) and the first identity in ( |3.7|) shows that these expressions agree. 

A similar procedure can be adopted to show that the formula in Q4.3j ) for S^{x,y) in the 
Laguerre case agrees with that implied by ( 4.26| ). Specifically, the task is to show 



S±{x,y) = §i(x,y) 



ca— >a— 1 



(A. 



This we do by comparing coefficients of y a / 2_1 e-»/ 2 L^- 1 (y), k = 0, . . . ,2N. In the case k = 2N 
on the l.h.s. use of (|A.5| ) shows this coefficient is equal to 



-x a l 2 e~ x 



(2N)\ 



;/2 (27V-l)!)^ JV _ 1 (x) 

T(2N)T(a + 2N) [ } 4T(2N + a) 



while on the r.h.s. we read off from ( f4.26| ), ( 3.24| ) and ( |3.4| ) that the same coefficient is 



(2N)\ 



4T(a + 2JV - 1) 



t a / 2 - 1 e-*/ 2 Z,S^- : L 1 (t) dt. 



One now checks that both the expressions vanish as x — > (for the latter this requires (A. 



and then verifies the equality of the derivatives. Equating coefficients of y a ^ 2 1 e y ^ 2 Lf,{y) for 
k < 2N gives an identity equivalent to the first formula in (3.7), thus verifying their equality. 
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